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We present a numerical algorithm to evaluate the self-consistent field theory for melts composed of 
block copolymers with multiply-branched architecture. We present results for the case of branched 
copolymers with doubly-functional groups for multiple branching generations. We discuss the sta- 
bility of the cubic phase of spherical micelles, the A15 phase, as a consequence of tendency of the 
AB interfaces to conform to the polyhedral environment of the Voronoi cell of the micelle lattice. 



I. INTRODUCTION 



Block copolymers provide an ideal route to engineer- 
ing well-controlled structure on nanometer length scales 
[l 0]. Through control over the chemical architecture, 
these systems can be tuned to self-assemble into periodic 
structures of an astounding variety. A plethora of new 
phases and structures have been identified in dilute di- 
block systems 0, triblock systems 0, and confined 
diblocks 1^ . One might think there is hardly more to say 
about melts of the simplest of block copolymer architec- 
tures, the neat linear AB diblock copolymer. It is well 
known that these linear diblock copolymers display a host 
of ordered phases: spheres, cylinders, lamella and the bi- 
continuous gyroid However, we have argued 0, 
that the tendency to minimize the AB interfacial area 
should stabilize a new cubic phase with Pm3n symmetry, 
the A15 lattice. The subsequent synthesis and characteri- 
zation of PEO-docosyl dendrimeric diblocks corroborated 
our prediction ITol | and was in agreement with the self- 
consistent field theory (SGFT) phase diagram for mik- 
toarm star copolymers 0. In this article, we provide 
the details of SGFT for branched architectures and, to 
our knowledge, the first SGFT phase diagrams for true, 
multiply-branched dendritic diblock copolymers. 

The serial development of new chemical synthesis 
routes is typically a costly and slow means for exploring 
the consequences of novel copolymer architectures. It is 
therefore desirable to develop theoretical tools for the effi- 
cient computation the phase behavior which can system- 
atically map out novel phase properties for a broad class 
copolymer architectures. Milner and Olmsted developed 
a strong-segregation theory (SST) approach to the phase 
behavior of A„B„j miktoarm star copolymer melts, appli- 
cable in the ~^ cxd limit, where x is the Flory-Huggins 
immiscibility parameter and N is the total number of 
chemical segments in the copolymer 0, ^| . For asym- 
metric copolymers, say for n > m, the effective spring 
constant of the more abundant polymer block is r? jrv? 
times more stiff than other block. Such asymmetry leads 
an enhanced stability of phases with strong interfacial 
curvature, and thus, spherical and cylindrical micelles 
are predicted to dominate the phase behavior for large 
molecular asymmetry 11] . Fredrickson and Frischknecht 
introduced an approximate SST approach to multiply- 
branched dendritic copolymers [13, and Pickett devel- 
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FIG. 1: A schematic of the branched molecular architecture. 
The first generation A-block contains fN segments. Each 
higher generation B-block is composed of (1 — f)N /Ms seg- 
ments. Here, the branching of each generation is 2. In the 
mean-field approximation, it is necessary to define only a sin- 



gle coordinate, So 
generation. 



for the set of branching points of the a 



oped a more refined self-consistent brush analysis for 
dendritic copolymer melts Both works showed a 

similar increase in stability of high interface curvature 
phases. Despite the analytic transparency of these SST 
calculations, the results of these calculations are predi- 
cated on many assumptions about the detailed structure 
of the micellar aggregates. In particular, certain assump- 
tions must be made concerning the interfacial shape and 
the direction in which copolymer chains stretch in the 
aggregates 0|. Because the free energy differences 
between phases are small, the presence of these undeter- 
mined degrees of freedom makes the task of locating the 
true free energy ground state analytically cumbersome, 
if not impossible. 

In Section II we present the theoretical derivation of 
the SGFT for multiply-branched copolymer melts from 
the full classical partition function of this system. We 
present an algorithm for the SGFT of block copolymers 
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within a specific class of multiply-branched architectures 
(see Figure nj. Like the SOFT approaches of Matsen 
and Schick for linear AB diblock copolymers and for 
(AB)„ starblock copolymers this approach makes 
no approximation beyond the approximation of mean- 
field in the monomer concentration profile. Therefore, 
this SCFT fully captures the copolymer chain fluctua- 
tions in the presence of the average concentration profile 
of constituent monomers. Moreover, this approach effi- 
ciently minimizes over all possible copolymer configura- 
tions, precluding the variational assumptions often neces- 
sary in the SST calculations. Finally, a numerical imple- 
mentation of SCFT is not limited to the infinite P^- 
rameter range. Given an arbitrary amount of computing 
time the equilibrium phase can be determined for any fi- 
nite value of xN. Practically, SCFT provides an efficient 
means of computing the mean field free energy of most 
phases for < 100 Tt^. In Section III we present the 
results of an application of SCFT to a series of branched 
copolymer melts within a specific class of this structural 
motif, specifically, copolymers which branch doubly with 
each successive generation. We discuss these results in 
context of elastically asymmetric copolymer melts and 
the geometry of the AB interface. We conclude in Sec- 
tion IV. 



II. SELF-CONSISTENT FIELD THEORY FOR 
BRANCHED ARCHITECTURES 

Our approach to multiply-branched diblocks is an ex- 
tension of the SCFT approach to linear diblocks and 
starblock copolymers pioneered by Matsen and Schick 
[TEl IT^ . While the derivation of the mean- field free en- 
ergy for the mulitply-branched system follows directly 
from the results for the linear and starblock architec- 
tures, we present its full derivation here since the sub- 
sequent evaluation of that free energy requires a slightly 
more generalized approach. Nevertheless, where possible, 
we attempt to keep the notation consistent with theirs. 

We consider a system of total volume, V, containing n 
branched copolymers. These copolymers are each com- 
posed of N total segments. Without loss of generality, 
we define the segment volume for both monomer type 
to be , so that the total volume of the system is 
V — — . The statistical segment lengths for the A and 
B-type monomers are denoted by a a and ob . The volume 
fraction of A-type monomer in the system is denoted by 
/. Thus, each chain is composed of fN A-type segments 
and (1 — f)N B-type segments. The architecture of our 
molecule is shown in Figure ^ The first generation is 
a single A-block. Grafted onto this are {g — 1) genera- 
tions of equal length B-blocks. The branching of the a'^ 
generation is given by rja so that the total number of B 
blocks, A/b, is given by, 

A^B =7/2(1 + ??3(l + f?4(... (! + %-! (1 + %))...) . (1) 

We define a coordinate along the polymer, s, so that 



within any chain portion of length As there are {As)N 
segments. Thus, in these coordinates, the length of the 
A block is given by Asa = / and that of the B sections 
is given by Asb = (1 - /)/AAb- 

A particular melt configuration is specified by n 
branched curves in space, r/3(s), the course-grained po- 
sition of the (sA^)**^ segment of the chain. At this 
point, we do not introduce an explicit parameterization 
of the full branched configuration. It suffices to demand 
that the first generation curve is joined to 772 second gen- 
eration curves which are each joined to 773 third curves, 
etc.. Given this set of branched curves, we define the 
dimensionless segment density operators. 



Po 



dsjis) Siv-vp{s)) , (2) 



{r)^-J2 ds (l-7(s)) S{r- 

po fr[ Jo 



, (3) 



where j{s) is a function which is equal to 1 when s lies 
along an A portion of the chain and when s is along a 
B portion of the chain, and the integration range is over 
the entire branched curve. In a neat system, the allowed 
melt configurations arc incompressible, and thus we are 
constrained to consider configurations for which 
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The full partition function for the melt is the functional 
integral over n branched curves: 
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where a normalization factor is absorbed into the func- 
tional measure, [dr^], r(s) = dr{s)/ds, k = ua/cib mea- 
sures the relative length of the A and B segments, and 
a = GA. The Flory-Huggins parameter, characterizes 
the repulsive interaction between unlike monomers. 

We can use the identity /[c?<i>A,B] (5[$A,B(r) — 
0A,B(r)] = 1 to transform into a functional integral 
over the monomer distributions. Introducing fields conju- 
gate to the total and individual segment concentrations, 
we have explicit representations of the delta- functionals, 

<5[l-(^A(r)-0B(r)] = 



m exp ^ / dr S(r)[l - Mr) - M^)] } , (6) 
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and, 

<5[«'a,b(i') - 0a,b(i")] = 

[dWA^B] expJ^ /"driyA,B(r)[$A,B(r)-0A3(i-) 



(7) 



where the hmits of integration of the conjugate fields are 
±ioo. Inserting these representations and the above iden- 
tity into ijnj and integrating over the delta functions in 
(O and (|3J), the full partition function is given by, 



Z = 
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where Q[Wa{^), VFB(r)] is the partition function for a sin- 
gle non-interacting, branched chain subject to the spatial 
field, WA(r) acting on first generation of the chain and 
WB(r) acting on the higher generations: 
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In general, it is not possible to evaluate the functional 
integrals in 0. Nevertheless, in the limit where N is 
large, fluctuation contributions to the partition function 
are small, and the integral is dominated by its saddle- 
point, where the free energy per chain, —^^^-\nZ, is 
minimal |l8l [T9| . The saddle-point approximation, of 
course, yields the mean-field results. 

To obtain the mean- field result, we solve for the field 
configurations, [0A(r), </>B(r), WA(r), WB(r), C(i")]7which 
minimize the free energy (that is, the lower-case func- 
tions are the extremal values of the upper-case functions). 
Minimizing with respect to <I'a(i"), $B(r) and S(r), re- 
spectively, we obtain the mean-field equations: 



wa{v) = x^0B(r) 



«;B(r) =x^0A(r)+C(r) 



<^A(r) + <^B(r) 



(10) 
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Minimizing with respect to WA(r) and WB(r), respec- 
tively, we find expressions for the mean-field densities. 
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where we have defined Q = Q[wA(r), WB(r)]. 

Upon inspection, it is clear how these relations consti- 
tute the mean-field theory result of the full problem. We 
have replaced the problem of multiply-branched chains 
mutually interacting, with the problem of non-interacting 
chains subject to the fields WA(r) and WB(r). These fields 
are chosen to represent the mean-field interactions pro- 
duced by the monomer distributions (/'A(r) and (/)B(r), 
That is, from ifTUI) and ((TT|l it is clear that A- type (B- 
type) monomers experience a repulsion proportional to 
XN times the local density of B-type (A-type) monomers 
and a repulsion due to the overall incompressibility of the 
system, given by ^(r). Because the mean- field incom- 
pressibility constraint, H12|l . depends only on the total 
monomer density, ^ (r) contributes equally to both poten- 
tials, WAiic) and WB(r). Hence, we see that ^(r) is simply 
the Lagrange-multiplier field which allows us to fix the 
combined, local segment concentration to po- Moreover, 
the average segment distributions, H13|l . are simply the 
average distributions produced by non-interacting chains 
subject to the fields u'A(r) and WB(r). Thus, eqns. itTUI) - 
H13|l provide a fully self-consistent set of equations, which 
can be solved to yield the mean-field result. Once the 
wa(j) and WB(r) are found, we can compute the mean- 
field free energy per chain, 

= - \nQ-V-^ [ dr[u.A(r)0A(r)+u;B(r)(^B(r)] 
+ [ dv xNM'^)^b{v) . (14) 



The first line of (|14|l gives the entropy per branched 
chain, and the second line gives the enthalpic, or inter- 
action, contribution to the free energy. 

For a given set of monomer potentials, WA(r) and 
WB(r), Q can be evaluated. Wc start by defining the 
Green's function, or propagator, for a continuous, un- 
branched portion of the chain. 
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where this path integral is carried out over all paths, 
r/3(s), which such that r^(si) — and Yp{sf) — r/. We 
absorb a normalization into [dTp] so that the integral of 
the propagator over the coordinates and is indepen- 
dent of arc length, Sf — Si. This is the same as demanding 
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FIG. 2: A schematic representation of the probability cap- 
tured by the end-distribution functions, g^(r, s) and q{r,s), 
for a 4-generation molecule. For the point, s, q\r,s) is pro- 
portional to the probability that the dashed portion of the 
chain has diffused to the position, r. For the same point, 
q{r, s) is proportional to the probability that the dotted por- 
tion of the chain has diffused to the the same position. The 
probability that the point is at r at s is the product of q and 



that the probability of any portion of this chain having 
any configuration (in the absence of external potentials) 
is independent of the number of segments it contains. 
Note that G{ri,Si;rf,Sf) is identical to the imaginary- 
time quantum mechanical amplitude (with s —^ —it) for a 
particle of mass, Na'^/3 (or Nq^/Sk^ when 7(5) ~ 0), in 
the potential, —wa{^) (or — ^^(r) when 7(5) = 0) mov- 
ing from Ti at the initial "time," Si, to at a later "time," 
Sf. Therefore, we know that G(ri, Si;r f, Sf) obeys the 
imaginary-time Schrodinger equation, or diffusion equa- 
tion and, unlike its interpretation in quantum mechanics, 
represents a probability and not an amplitude. We make 
explicit use of this fact below. 

To capture the branched architecture of the chain we 
define the end-distribution functions. These functions 
compute the statistical weight of a chain diffusing along 
its trajectory to some position in space. That is, we 
define a function, (r, s) , which is proportional to the 
probability that the branched chain diffused from one 
of its free ends at Sg, where Sa is the length coordinate 
corresponding to the branching point of a^^ generation 
(see Fig.mi. Note that for Sg_i < s < Sg this function is 
simply the probability that an unbranched chain diffused 
from its free end to s at some position, r. But if Sg-i < 
s < Sg-2! then q'^{r,s) is proportional to the probability 
that r]g free ends diffused from Sg to some intermediate 
position, say rg_i, at Sg_i and then diffused on to r at 
s (see Fig. |2l. Thus, as s decreases towards the free 
end of the A block at sq, ^^(r, s) assumes the probability 
of all the higher generations diffusing "into" this lower 
generation branch. We will refer to this diffusion from Sg 
towards sq as "backward motion." Note that in terms of 



the probability distributions all paths diffusing from any 
of the rjg free ends are equivalent, and therefore, q^{r, s) 
is well defined. 

We summarize the above definition by writing q^{r, s) 
in terms of our unbranched propagator, G(ri, s^; r/, s/): 



/ dva, G(r, 



for Sq-I > s > So 



(16) 



where q^{ra, ) indicates that we take the value of this 
function from the end of the higher generation at (just 
before the branch point). If we normalize our propagator 
so that hm^^^si G(r,i, Si; r/, s/) — S{rf — r^), then we 
establish a set of boundary conditions for q^{r,s) at its 
free end at Sg and each branching point. 
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<zt(r,,s-) = [gt(r„,4)] 
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where (/^(r, s~) is the limit of the function as s ap- 
proaches Sa from below (just after the branching point). 
Thus, at a given branching point, Sq the value of ^^(r, s) 
changes discontinuously, from q^r, s+) to q^{r, s~), since 
the function assumes the probability of the other higher 
generation branches meeting it at that point. 

Because q^ (jc, s) is defined in terms of the propagator, 
G{ri, Si;rf, Sf), we know that it will obey the same dif- 
fusion equation as the propagator. Namely, 



^^V^gt _ ii;^(r)(jt ^ for So < s < si , 



ds 



V^gt — WB{r)q^ , for si < s < Sg 



(19) 

It should be understood that we will solve these first- 
order equations for the unbranched portions of the chain 
and use the branching points to determine boundary con- 
ditions; hence, we do not need to worry about differenti- 
ating at branching points. 

Because Eq. p9f) is a linear equation for q^ which is 
first order in s, given any set of fields, ^^(r) and ^^(r), 
we can solve for q^{r,s) for all segments. First, using 
(|19|l and H17|) we solve for the q'''{r,s) for the g^^ gener- 
ation. Then, we can use H19|l . (|18|l and our solution for 
q^r, Sg-i) to solve for the {g—1)^^ generation of g^(r, s). 
Likewise, we can then iteratively solve for all lower gen- 
erations until we get to the first. 

Once the value of q'^{r, s) is known for all s down to sq, 
we can compute the single-chain partition by integrating 
this backward motion end-distribution function over the 
position of the free end of the A-block, 



Q= dv gt(r,so) 



(20) 



However, in order to compute the mean-field melt free 
energy we need to calculate the average monomer dis- 
tributions, 0A(r) and 0s (r), created by the monomer 
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potentials, wa{y) and wb{j^)- By introducing another 
end-distribution function, g(r, s), we can compute the 
functional derivative of — In Q with respect to these fields 
directly. 

We define q{r, s) to be proportional to probability that 
a chain configurations diffuses in the "forward" direction 
from its other free end (the free end of the A block at sq) 
along one of the branched trajectories of the molecule to 
s at the position r (see Fig. IJJ . At branching points, Sa , 
q{r, s) assumes the probability that {rja+i — 1) branches 
have also diffused from their free ends at Sg to Tq, at Sq. 
This is to say that q{ra, sj^) contains not only the prob- 
ability that the sq end diffused to this point but also the 
probability that all of the other branches, not including 
the currently diffusing path, have diffused to to Tq at 
to meet it. This property makes q{r, s) convenient 
for computing the average monomer distributions. Using 
the above definition we have, 

„r..,=/*„G,r„..,;r,,),(r„,;,|,.(r<.,4)r- 

for Sa > S > Sa + 1 . (21) 

The corresponding boundary conditions for q{r, s) are 
given by, 

q{r,s+)^l, (22) 



g(r,4) = g(r,s-) [gt(r„, s+)p-i . (23) 

Since the "motion" of the diffusion along the chain is 
reversed from that of q'^{r,s) the diffusion equation for 
<;(r, s) is the same as (IT^ except with a plus sign ap- 
pearing on the left hand side. In analogy with gt(r, s), 
we must first solve the diffusion and (|22|l for the first 
generation of g(r, s). We then use our second generation 
solution of g'l'(].^g) and the first generation solution of 
q{r, s) in to find the solution for the second genera- 
tion. We can repeat the process to solve for q{r, s) over 
the entire length from sq to Sg. 

It is not difficult to show that the monomer distribu- 
tions, given by Eq. H13() . can be computed by, 

0A(r) = -^ / 'rf^9(i',s)9^(r,s) , (24) 

J So 



0B(r) = -^Yl -^B.a / ds q{r, s)q^r, s) , (25) 
where V = — and Afn a is the number of B-blocks in the 

Po ' 

Q,th generation, which is simply given by rjaTja-i . • .?72. 
Thus, the mean-field free energy, ((Tl|l . can be computed 
entirely with the end-distribution functions, ^(r, s) and 
g'l'(r, s) . 

While real-space methods for numerically solving these 
diffusion equations exist, [T^ l20j | these methods tend to 



be computationally intensive for melt phases with spatial 
variation in three dimensions. Instead, we use Fourier ex- 
pansions of the functions to solve for q^(r, s) and (/(r, s) 
given an arbitrary set of external fields, u;^(r) and ^^(r). 
Since we know that equilibrium structures are themselves 
infinitely periodic structures, we expect that we can very 
accurately describe mean-field results with a finite num- 
ber of Fourier terms included the expansion. For up to 
moderately large degrees of segregation (for ^ 50) 
the spectral methods of |0| and 0| allow for the rapid 
and very accurate exploration of mean-field thermody- 
namics T^ . We present the spectral form of our SCFT 
for multiply-branched copolymer melts in the Appendix. 

III. DOUBLY FUNCTIONAL BRANCHING: 
THE ROLE OF INTERFACES 

Using the SCFT derived in the previous section we 
computed the ^ 40 mean-field phase behavior for 
multiply-branched copolymer melts where the branching, 
or functionality, of each generation is 2. We compute 
the phase behavior for g = 2 . . . 6 for monomers of equal 
segment size, k = 1. To achieve a precision of ±10"'^ in 
/ and ±10~^ in x^ require a precision in the free 
energy of better than ±10"'*. This requires the use of up 
to 908 basis functions for some phases. The mean-field 
phase diagrams for 3 < g < 6 are shown in Figs. Eland 
0] We have already reported on the phase behavior for 
5 = 2, the AB2 miktoarm star 

The thermodynamics of these melts are strongly influ- 
enced by the introduction of the multiply-branched ar- 
chitecture. Compared to the predicted phase behavior of 
linear AB block copolymer melts, the phase boundaries 
of these branched copolymer melts are skewed system- 
atically towards larger values of / for most phases [isj. 
This indicates an enhanced preference for phases where 
the branched polymer domain is on the convex side of 
curved AB interfaces. In Figure we plot the strong- 
segregation [xN = 40) phase boundaries as a function of 
branching generation. The preference for morphologies 
with the branched, B domain on the outside of a highly 
curved interface is increases with increasing generation. 
For example, spherical micelles where the A blocks form 
the core region are stable up to / = 0.275 for g — 2 but 
stable up to / = 0.350 for g = 6. This effect is well 
established for copolymer architectures with elastically 
asymmetric blocks [Tl Illl2lll23| . 

In general, elastic asymmetry stems chiefly from two 
sources — asymmetric monomer sizes and asymmetric 
copolymer architecture. Milner demonstrated within 
SST that the elastic asymmetry between copolymer 
blocks of a A„B„i miktoarm star is captured by the pa- 
rameter, e — ^( y^'^: where and are the re- 
spective volumes of the A and B segments . From this 
analysis it can be shown that the effective spring constant 
of the B brush domain is a factor of times the value of 
the symmetric case (for e = 1). For e > 1, the molecu- 
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FIG. 3: Phase diagrams for g = 3 and g = 4. Dis labels re- 
gions where the melt is disordered. Stable regions of ordered 
phases are labeled: (Lam) lamellar; (Gyr) gyroid, laSd sym- 
metry; (Hex) hexagonal-columnar, pQmm symmetry; (A15) 
sphere phase, PmSri symmetry; (BCC) body-center cubic lat- 
tice of spheres, Im3m symmetry; and (FCC) face-centered 
cubic lattice of spheres, Fm3m symmetry (2Jj. The circle 
marks the mean field critical point through which the system 
can transition from the disordered state to the Lam phase via 
a continuous, second-order phase transition. All other phase 
transitions are first-order 



lar asymmetry leads to the stabilization of morphologies 
where the B polymer block composes the outer corona 
of spherical and cylindrical domains for larger values of 
A composition than is observed for elastically symmetric 
copolymers psj . 

It is desirable to have a similar quantitative mea- 
sure of the elastic asymmetry for copolymers with this 
multiply-branched structural motif. However, in con- 
trast to the miktoarm star architecture, the elastic en- 
hancement of multiply-branched domains depends on the 
aggregate morphology. Using the Alexander-de Gennes, 



FIG. 4: Phase diagrams for g = 5 and g — 6. Labels appear 
as in Figure 121 

strong-segregation analysis employed by Frischknecht 
and Fredrickson we find, for example, that the stiffness of 
a lamellar B domain in these doubly-functional copoly- 
mer melts is enhanced by a factor of 4(8^"^ — 1)/[7(2S^^ — 
1)] over the linear, unbranched case This corre- 

sponds to factors of 4, 12, 2|2 ^ 41,7^ 155 and 604 multi- 
plying the stretching free energy of a lamellar B domain 
for the g = 2,3,4,5 and 6 case respectively. Pickett 
demonstrates, however, that the Alexander-de Gennes 
approximation provides an overestimate of the branched 
chain free energy whose error grows quickly with the 
branching generation . Based on the analysis of Pick- 
ett for a slightly different copolymer architecture we 
might expect that by relaxing the constraint that the 
chain ends are held at the tips of the brush, the SST 
stretching free energy of the branched B domain can by 
relaxed from the Alexander-de Gennes upper limit by 
factors of roughly 2.6, 5.5, 11.7 and 24.6 for g = 3,4,5 
and 6, respectively (the g = 2 case corresponds the AB2 
miktoarm star). This allows us to estimate more realis- 
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FIG. 5: The SCFT phase boundaries computed at xN = 40 
for 2 < g < 6 are depicted as open circles. For comparison the 
xN = 40 phase boundaries for linear diblocks are indicated of 
the / axis. Note the absence of a stable A15 phase for linear 
AB diblock copolymer melts. 



tic values of the elastic asymmetry in the lamellar mor- 
phology: 4 for g = 2; 4.6 for g = 3; 7.7 for g = A; 
13.4 for g — 5; and 24.6 for g — 6. While these are 
somewhat crude estimates, they provide reasonable cor- 
respondence between the SST phase behavior of these 
multiply-branched copolymers and AB„ copolymers |^. 
For example, the respective xiV = 40, Gyr-Lam and A15- 
Hex transitions occur at / = 0.546 and / = 0.349 for a 
melt of AB5 miktoarm stars, corresponding to an elastic 
asymmetry of 25. This should be compared to the same 
transitions which occur at / = 0.550 and / — 0.350, 
respectively, for a 17 = 6 branched copolymer, with an 
estimated elastic asymmetry of 24.6. 

We note the stability of the cubic A15 phase in these 
melts. We have argued 3 that as ~^ 00 and in 
the limit that the AB interface of a sphere phase is con- 
strained to adopt the same shape as the lattice unit cell 
that the A15 should be the equilibrium structure. In this 
limit the relative stability of competing arrangements of 
micelles can be assessed purely in terms of two geometric 
moments of the Voronoi polyhedra of the lattice: the area 
and the second- moment of the lattice . If Rx{^) mea- 
sures the radial distance from the center to the surface 
of the Voronoi cell of lattice X at solid angle, tt, then we 
can compute the area in terms of the area of a spherical 
cell of equal volume, 



A 



X 



1 jdn[Ri{Ci) + iv^Rxmy 



(47r)i/3 



2/3 



(26) 



where — 9-^ -f We can also define the second- 



moment, or "stretching" moment, of the Voronoi cell. 



Ix = (4^) 



2/3_ 



/ dhR\{n) 
{jdnRj.mf 



(27) 



where again we have normalized by the same measure 
for a spherical cell of equal volume. It can be shown that 
the free energy per chain in a micellar phase arranged 
in lattice X is simply given by Fx — Fo(yi^X)^/'^, where 
Fq is the free energy per chain for the case when the 
Voronoi cell is approximated as a sphere 

Given 

these geometric measures for all candidate arrangements 
of spherical micelles we can assess the relative stability 
of these phases in this limit where AB interface has the 
same shape as the unit cell of the lattice. It was discov- 
ered by Weaire and Phelan that the space partition of the 
A15 lattice has the lowest area of all known equal volume 
periodic partitions of three dimensional space ^231. It is 
for this reason, despite the fact that the BCC lattice has 
a smaller second moment, that the A15 lattice is most 
stable among the lattice arrangements of spherical mi- 
celles when AB interfaces have adopted the shape of the 
Voronoi cell in which they are confined. In particular, 
this limit predicts that the free energy per chain for the 
A15 phase is 0.14% and 0.61% lower than the BCC and 
FCC phases, respectively. Of course, there are finite x^ 
corrections to this asymptotic limit due to chain fluctua- 
tions which are neglected in the strong-segregation limit, 
but the lowest order corrections which distinguish be- 
tween morphologies are smaller than the leading order 
free energy term by a factor of (xiV)~^/^ '26!, '27']. 

The conclusions of our SST analysis are valid in the 
limit that the AB interface has adopted the polyhe- 
dral shape of the lattice Voronoi cell. It is well known 
that constraining a micelle to occupy a polyhedral unit 
cell frustrates the internal configuration of the aggregate 
[28L l2^ . Chains which extend along directions towards 
corners of the Voronoi cell must stretch further than those 
extending towards the walls. The difference in tension in 
these chains leads to a tendency to distort the AB in- 
terface from its ideal, uniformly curved shape into the 
polyhedral shape of the Voronoi cell. Of course, the mi- 
celle will adopt some compromise between the uniform 
curvature and relaxed outer chain stretching which will 
be determined by the relative importance of outer chain 
stretching and the forces which pull inward on the AB 
interface, namely the surface tension and the inner chain 
stretching. We demonstrated how the tendency for 
cylindrical micelles in the Hex phase to adopt an hexag- 
onal interface shape is enhanced both by an increase in 
inner domain volume fraction, /, and the elastic asymme- 
try between the coronal and core polymer domains, e. In 
particular, wc found that although AB interfaces in mi- 
celles for symmetric molecules (e.g. linear AB diblocks) 
remain relatively unperturbed by the lattice symmetry, 
cylindrical micelles composed of very elastically asym- 
metric copolymers (e ^ 3) have interfaces which are very 
nearly hexagonal in regions where the Hex phase is ther- 
modynamically stable. Given the stability of the A15 
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FIG. 6: A view of the Voronoi cell of the BCC lattice, a 
truncated octahedron, with half of one hexagonal face and a 
quarter of one square face removed to reveal the inside. The 
edges are drawn as black lines, the outside is shown as blue 
and the inside is shown as yellow. The vectors connecting the 
center of the cell to the corners along the (210) directions and 
the nearest walls along the (111) directions are shown. 
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FIG. 7: Plot of the measured distortion, a as defined in eqs. 
ll^ and l^nj, measured from SOFT results for the BCC phase 
of branched copolymer melts for generations, 2 < g < 6. 
For comparison the dashed line shows the same distortion for 
linear AB diblock copolymer melts. 



phase in the present system (Figure EJ, it must be that 
the interfaces of the sphere phases are also substantially 
distorted by the polyhedral environment of the lattice 
Voronoi cell. 

We can quantify the extent to which AB interfaces in 
sphere phases adopt the shape of the Voronoi cells from 
our SCFT results. A measure of the distortion of the 
interface of a micelle in the BCC phase from the ideal 
spherical shape is the difference of the distances from 
the center of the micelle to the interface along directions 
towards the closest face of the Voronoi cell, the (111) 
direction, and towards the corner of the Voronoi cell, the 
(210) direction, 

^^^(210)^^^ (28) 
R{2w) + -R(iii) 

where i?(2io) and i?(iii) are the radial distances to the 
AB interface along those directions (see Figure IHI). For 
a spherical interface we have Ssph = and for an inter- 
face which as the truncated-octahedron shape of the BCC 
Voronoi cell, Sbcc = (%/5 - \/3) / ( %/5 + %/3 ) ~ 0.127. By 
normalizing measured values of S by Sbcc, we can assess 
the polyhedral distortion on the scale set by the shape of 
the BCC Voronoi cell. Therefore, we use, 

a^j^ (29) 
obcc 

to quantify the polyhedral distortion of the interface as 
a function of molecular architecture. Figure plots the 
shape parameter, a, measured from SFCT calculations 
for the BCC phase as function of / and branching genera- 
tion g. It is clear that the packing frustration introduced 
by the polyhedral Voronoi cell increases as the volume 
fraction of the core of the micelle grows. Although the 



close-packing limit of hard spheres in a BCC lattice is at 
a volume fraction of \/3n/8 ~ 0.68024, the cores of the 
micelles are highly deformed for / well below this. Since 
the outer chain stretching is responsible for this poly- 
hedral distortion, the tendency to adopt the truncated- 
octahedral shape of the BCC lattice is enhanced as the 
stiffness of the coronal region is increased by molecular 
branching. 

To further visualize this distortion we compute the 
mean-curvature, H, of AB interfaces extracted from 
SCFT results for melts at = 40 at the phases bound- 
ary between spherical and cylindrical phasess, the A15- 
Hex boundary, for g — 6, g — 2, and for linear diblocks. 
The distortion in these interfaces corresponds to mea- 
sured values of a = 0.32, a = 0.128 and a = 0.011, re- 
spectively. These surfaces are displayed in Figure |S1 and 
following the analysis of Matsen and Bates the inter- 
faces are shaded according to the local mean-curvature. 
As the polyhedral distortion increases, the deviation from 
constant mean curvature grows. The surface tension, 7, 
associated with an interface between unlike polymer melt 
domain scales as x^/^ 113. A patch of area, dA, of a 
curved interface experiences a force due the surface ten- 
sion which is given by 2H^dA j^. Since these micelle 
configurations are saddle-points of the free energy, we 
know that the force due to the tension pulling inward on 
the interface is balanced by a net force pulling outward 
on the micelle interface, that is, the interface must be in 
mechanical equilibrium. In copolymer micelles the com- 
pensating forces are due to a difference in the tension of 
the chains in the core and coronal domains. Therefore, 
variation of the mean curvature of the interface provides 
a direct measure of the chain tension pulling on the inter- 
face. Regions of high interfacial curvature, towards the 
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FIG. 8: The AB interfaces extracted from the SOFT calculation (the surfaces at which 0a (r) = 0.5) for the BCC phase along 
the thermodynamic phase boundary separating spherical and cylindrical morphologies: (a) for linear diblocks at / = 0.166; (b) 
for g = 2 branched copolymers at / = 0.275; and (c) for g = 6 branched copolymers at / = 0.350. The interfacial distortion 
at these points corresponds to measured values of a = 0.011, 0.124 and 0.321, respectively. The surfaces are shaded according 
to the local mean curvature, H, measured in units of the average mean-curvature, (H). The variation of the mean-curvature 
provides a direct measure of the variation of the polymer chain tension at the interface, due to the polyhedral environment of 
lattice Voronoi cell. The standard deviation of the curvature, aH, for each surface is given in units of (H). 



edges and vertices of the BCC Voronoi polyhedron, cor- 
respond to regions where the coronal B domains pull rel- 
atively strongly on the interface. Conversely, relatively 
flat regions on the AB interface, towards the nearest- 
neighbor faces of the polyhedron, indicate that the chain 
stretching is relatively low. 

IV. CONCLUSION 

From our analysis we see that for elastically asym- 
metric copolymers the polyhedral shape of the lattice 
Voronoi cell forces the micelle configuration to deviate 
drastically from the limit of uniform interfacial curva- 
ture. While we have argued that in the limit of perfectly 
polyhedral interfaces the A15 phase is most stable among 
the sphere phases and that very elastically asymmetric 
micelles approach this polyhedral limit in regions where 
sphere phases are thermodynamically stable, it has yet 
to be shown that for small distortions (for a < 0.35) the 
minimal Voronoi cell area argument should apply. For 
example, the reduced area of the g = 6, BCC interface 
of Figure |H1 is 1.0094, to be compared to the reduced 
area of the truncated-octahedron of the BCC Voronoi 
cell, 1.0990. In this sense, AB interfaces of physical mi- 
celles seems to be distorted less than about 10% towards 
their Voronoi polyhedra. Nevertheless, we argue that 
the polyhedral interface limit of the micelle configura- 
tions sets the scale of the frustration. While the true 
micelle interfaces are some interpolation between a spher- 
ical and polyhedral shape, the scale of the frustrated free 
energy is set by the polyhedral interface upper bound. 
As mentioned above, in the limit of polyhedral interfaces 



as xiV CX3, SST predicts that Fbcc = 1.0014^^15 and 
Fpcc — 1.0061i^Ai5. We can compare this prediction to 
the results of our SCFT calculations along the Hex-Al5 
phase boundary at = 40 for g = 5 branched copoly- 
mers (at / = 0.340) for which we find Fai5 = 6.296nfcsT, 
Fbcc = 6.3UnkBT and Fpcc = 6.326nkBT, corre- 
sponding to 0.28% and 0.46% higher free energy than the 
A15 phase for BCC and FCC phases, respectively. On 
the scale of these small free energy differences, the anal- 
ysis of our geometrical limit is a necessary component 
of any rational explanation for lattice choice. Therefore, 
while such a calculation has yet be carried out, we expect 
that a more detailed SST analysis of the relaxed configu- 
rations of micelle interfaces will bear this argument out. 
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APPENDIX: SPECTRAL SCFT 

Following Matsen and Schick 0, we define a set 
of orthogonal basis functions, fi{r), which have the pe- 
riodic symmetry of our copolymer phase. We expand all 
of the necessary functions of position in this basis, so 
that g{r) = J2i 9ifii^)- For example, the BCC phase of 
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spheres can be described by the set of functions with 
/m3m symmetry. The functions are normahzed such 
that, 



dv f,{r)fj{r) = % 



(A.l) 



In addition, we demand that these functions are eigen- 
functions of the Laplacian operator so that, 



(A.2) 



where D is the length scale of the periodicity of the sys- 
tem. The set of functions is ordered in an increasing 
sequence in Ai, and Ai is set to (or /i(r) = 1). Be- 
cause the product of two basis functions, /i(r)/j(r), has 
the all symmetries of the basis, it belongs in the same 
space of functions as our basis. Thus, we can write the 
product as expansion in our basis functions. We define 
the coefficients, Tijk, of the expansion so that. 



/.(r)/,(r)=^r,,,./4r) . 



(A.3) 



Alternately, given the set of basis functions invari- 
ant under all of the symmetry operations of the 
group, this coefhcient can be computed by Tijk = 

With these definitions and a finite Fourier expansion 
of all functions of position we can rewrite the diffusion 
equation, (17), as a matrix equation. 



r Ej ^y ?]. for So < s < si , 
I J2j Bijq], for si< s < Sg , 
where we have defined the matrices, Aij and Bij, 



(A.4) 



ijk 



k 



(A.5) 



(A.6) 



Since these are symmetric, real matrices we can 
diagonalize Aij and Bij by orthogonal transforma- 
tions, such that, Y.k.i\0'A\ikAki[OA]ij = AiSij and 
J2k.i\-^B^ikBki[OB]ij = B^Sij, where Ai and Bi are the 
eigenvalues of Aij and Bij, and [OA]ij and [Ob]^ are 
the orthogonal matrices which diagonalize, Aij and Bij, 
respectively. The matrix, 

tI..{s'^s) = ^[OA]a.exp{-A(s'-s)}[0^],fe , (A.7) 

k 

transfers the A-block solution to (|A.4|I from s to s', and 
the matrix Tg jj(s' — s) does the same for the B-block 



solution. Using these matrices we can write the solutions 
for g|(s), 

{Ej Tlijis - si)A](si) , for So < s < si , 
J2j TI ij(s - SQ)A](Sa) , for Sa-l < S < Sa , 
(« ^ 1) , 

(A.8) 

where Aj(sQ) are the boundary conditions for g|(s) at 
Alisc) = I dv [q\v, stW-^^h{v). (A.9) 



Thus, we have that Aj(sg) — Ja- 
in order to compute h\{sa) for the lower generations, 
we define the function, 

^^r\s») ^ I dv [<zt(r,4)]'"/,(r). (A.IO) 



Given this definition we have ?a|"'^''(sq) — (7j(s+) and, of 
course, l^\{sa) = ^'f'^^^\sa)- Using (|A.3p and the fact 
that the Fourier expansion of q\v, s+) = q\{s^)fi{v), 
it can be shown that, 

^^\^^) -Y.^^Jk q]{st) ^i"-^'(s„). (A.ll) 

In order to find q\{s) for the [g — Vf^ generation, we 
first compute q\{s'^_i) by HA.8|) and A|(sg) = Sn. Using 

(m) , 



(|A.lip we can then iteratively compute V'j (sg-i) for 
all m up to rjg^i. Then we will have Aj(sg_i) and q\{s) 
for Sg-2 < s < Sg-i. We can repeat this procedure until 
we have qj{s) down to the first generation. At this point 
we have computed the probability of a chain diffusing 
from its branched, B-block tips down to the end of A- 
block. That is to say, we have computed the single-chain 
partition function, Q/V ~ ql{sQ). 

To find qi{s) we have solve the same matrix equation 
as (|A.4|I except with a plus sign on the left-hand side. 
Again, the transfer matrix for the "reversed motion" A- 
block solution is defined by, 

TA,ij{s' - s) = J2[0AUeMMs' ~ s)}[OA]jk , (A.12) 



and TB,ij(s' — s) is defined similarly. The solution for 
qi{s) is 



qi{s) 



J2j 7A,ij(s - si)Aj(si) , for So < s < si , 

J2j ^B,ij(s - Sq)Aj (Sq) , for Sa-l < S < Sa , 

(a ^ 1) , 

(A.13) 
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where Ai{sa) are the branching point boundary condi- 
tions, 

A,(s„_i) = J dv q{r, s-_,)[gt(r, st^,)^--^ Mr). 

(A.14) 

Similar to the identity l|A.ll(l . it can be shown that 



(A.15) 



We use the fact that Ai(so) = ^n, the boundary condi- 
tion for the So free end, to compute qi{s) from l|A.13() for 
the first generation. We can then use (|A.15|I along with 

qi{si) and ipi^^~^\si) to find qi{s) for the second gen- 
eration. Repeating this process we can find Ai{sa) and 
and qi{s) for all generations. 

Using both qi{s) and q|(s) we can use and to 
compute the Fourier amplitudes of the monomer concen- 
trations, 4)A.i and 4>B,i, 



qliso) J so J 



ds ^gj(s)q^.(s)r 



ijk 



(A.16) 



^7VB,a / ds^qj{s)ql{s)rijk ■ 

(A.17) 



Finally, the single-chain partition function, Q, and the 
Fourier densities are used to compute the mean-field free 
energy per chain (up to an additive constant) are: 



= -\nq\{sa)-xNY,^A,^(^}B,^ ■ (A.18) 



ukb 



Now we need only to find the Fourier components of the 
field configuration, 'WA,i and WB,i^ so that we satisfy the 
self-consistency relations: 



WA,i - WB,i = xN{(l)B,i - (l>A.i) , 



Sil = 4>AA + ' 



(A.19) 



(A.20) 



Clearly, we have (f)A,i = / and <f)B.i = (1 — /), by def- 
inition. Moreover, we can set wa,i — X^{^ ~~ /) and 
wb.i — X^f- Since the monomer distributions, cjjA.i and 
(j)B,i, depend functionally on the fields, waa and wba, 
through Q[wA,i,WBA], Eqs. HA.19|I and (|A.20|I present a 
complicated set of non-linear equations. These are most 
easily solved by computing (j)AA and 4>B,i, for some ini- 
tial guess of WA,i and WB,i- Then waa and wba can 
be adjusted towards a solution of (jA.19|l and (|A.20p . 
The computation then proceeds iteratively until the self- 
consistent solution is found. 
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